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Abstract —In massive multiple-input multiple-output (MIMO) 
systems, it may not be power efficient to have a high-resolution 
analog-to-digital converter (ADC) for each antenna element. 
In this paper, a near maximum likelihood (nML) detector for 
uplink multiuser massive MIMO systems is proposed where each 
antenna is connected to a pair of one-bit ADCs, i.e., one for 
each real and imaginary component of the baseband signal. 
The exhaustive search over all the possible transmitted vectors 
required in the original maximum likelihood (ML) detection 
problem is relaxed to formulate an ML estimation problem. 
Then, the ML estimation problem is converted into a convex 
optimization problem which can be efficiently solved. Using the 
solution, the base station can perform simple symbol-by-symbol 
detection for the transmitted signals from multiple users. To fur¬ 
ther improve detection performance, we also develop a two-stage 
nML detector that exploits the structures of both the original ML 
and the proposed (one-stage) nML detectors. Numerical results 
show that the proposed nML detectors are efficient enough to 
simultaneously support multiple uplink users adopting higher- 
order constellations, e.g., 16 quadrature amplitude modulation. 
Since our detectors exploit the channel state information as part 
of the detection, an ML channel estimation technique with one- 
bit ADCs that shares the same structure with our proposed nML 
detector is also developed. The proposed detectors and channel 
estimator provide a complete low power solution for the uplink 
of a massive MIMO system. 


I. Introduction 

Massive multiple-input multiple-output (MIMO) is a trans¬ 
mission technique for cellular systems that leverages a large 
number of base station antennas to support many single 
antenna users IT]. With enough antennas, massive MIMO 
can eliminate inter-user interference completely using matched 
beamforming for downlink and matched combining for uplink 
if the base station has full channel state information (CSI). 
Because of its simple beamforming and combining structures, 
massive MIMO will be beneficial not only for cellular sys¬ 
tems but also for other wireless communication systems, e.g., 
vehicular-to-everything (V2X) communications using many 
antennas G). 

There are several practical constraints that are encountered 
when implementing massive MIMO systems. Because of the 
large number of antennas, it may not be possible to deploy 
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expensive and powerful hardware with small noise and dis¬ 
tortion at the base station. Prior work studied the impact of 
several hardware impairments including phase-drifts due to 
non-ideal oscillators and distortion noise caused by analog-to- 
digital convertors (ADCs) for massive MIMO systems. It was 
shown in El-El that having a large number of antennas helps 
to mitigate these hardware impairments, which confirms the 
benefit of massive MIMO. 

In this paper, we focus on uplink multiuser massive MIMO 
systems using extremely low-resolution of one-bit ADCs at 
the base station. Because the power consumption by ADCs 
grows exponentially with their resolution level 0, Q, using 
one-bit ADCs may be a practical way of implementing cost- 
efficient and green massive MIMO systems. It is expected 
that using one-bit ADCs is particulary beneficial for wideband 
communication systems (when properly handling frequency 
selectivity) that require high sampling frequency, which will 
become common in millimeter wave (mmWave) communica¬ 
tion systems il. Adopting one-bit ADCs for uplink massive 
MIMO, however, is challenging because of the severe thresh¬ 
old applied to the received signal. 

Using low-resolution ADCs for wireless communications 
has been investigated under various assumptions. It was shown 
in ||9l for one-bit ADCs and in ifTOl for low-resolution (one to 
three bits) ADCs that the capacity maximizing transmit signals 
for single-input single-output (SISO) channel are discrete, 
which is different from the unquantized output case. It was 
shown in CD that even a low-complexity suboptimal QPSK 
detector with one-bit ADC suffers only 1-3 dB signal-to-noise- 
ratio (SNR) loss compared to the unquantized case for SISO 
channels. In OH, El, it was shown that oversampling re¬ 
covers some of the loss in the SISO channel capacity incurred 
by one-bit quantization. The mutual information of the MIMO 
channel with quantized output was studied without optimizing 
input distribution in ifTD for one-bit ADCs and in lITSl . ifT^ for 
low-resolution ADCs. The input distributions were optimized 
to maximize the achievable rate of the quantized MIMO 
channel using one-bit ADCs in CD, El. Interestingly, the 
impact of having low-resolution ADCs is not that severe. For 
example, M showed that the rate loss due to using one-bit 
ADCs is 1.5793 bits/s for 4x4 MIMO with a quadrature phase 
shift keying (QPSK) constellation. It was proven in iflTl that 
the mutual information of MIMO with one-bit ADCs is only 
2 /tt times smaller compared to that of the unquantized MIMO 
case in the low SNR regime. With an optimized threshold, this 
gap can vanish in the low SNR regime EH. In Ha, digital and 
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analog combiners using low-resolution ADCs were compared 
in terms of the achievable rate for point-to-point mmWave 
systems where the derived achievable rates were tight only 
for the low SNR regime. 

The aforementioned work on quantized MIMO na-Eoi 
was restricted to point-to-point communications and focused 
on studying capacity or designing optimal transmit signals for 
quantized channel outputs using low-resolution ADCs without 
proposing specific signal processing algorithms for detecting 
the received signals. For point-to-point communications with 
low, but more than one-bit resolution (e.g., 2-3 bits) ADCs, a 
modified minimum mean square error (MMSE) detector was 
proposed in ED and later extended to an iterative decision 
feedback equalizer in ll22l . It is possible to extend the iterative 
detection technique developed in to multiuser scenarios. 
The iterative detector proposed in E2l exploits the assumption 
of nearest neighbor error of adopted constellation, which may 
not be effective for practical channel models. 

In this paper, we propose a practical near maximum like¬ 
lihood (nML) detector for uplink multiuser massive MIMO 
systems with one-bit ADCs. The complexity of the maximum 
likelihood (ML) detector grows exponentially with the number 
of users, which makes it difficult to use for massive MIMO 
with multiple users. The proposed nML detector, however, can 
be implemented with standard convex optimization techniques 
with marginal performance degradation when the number of 
receive antennas at the base station is large, and supports 
multiple users with arbitrary constellation sizes. To further 
improve performance, we also propose a two-stage nML 
detector, which exploits the structures of both the original ML 
detector and the one-stage nML detector. Numerical results 
show that, with marginal additional complexity, the two- 
stage nML detector can significantly improve the detection 
performance and achieve almost the same performance with 
the original ML detector. 

There are several recent papers on massive MIMO with 
low-resolution ADCs. It was shown in E3i and 12^ that 
linear detectors with one-bit ADCs work well for multiuser 
scenarios with QPSK and 16 quadrature amplitude modula¬ 
tion (QAM) constellations, respectively. An uplink multiuser 
massive MIMO detector was developed in ll25l but the detector 
was based on several bit ADCs and developed for the spatial 
modulation transmission technique ll2^ . Lor general symbol 
transmission in uplink massive MIMO systems with one-bit 
ADCs, a message-passing algorithm-based multiuser detector 
was proposed in 1271 for special constellation structures. The 
detector was later extended to arbitrary constellations in l28l . 
Our numerical studies show that our nML detector works well 
with both perfect CSI and channel direction information (CDI) 
while the detector in 1^ experiences the performance degra¬ 
dation with CDI when SNR is low. In 1^ . a multiuser detector 
using low-resolution ADCs was developed based on convex 
optimization, which has a similar structure with our nML 
detector. The convex optimization process in l29l . however, 
is optimized separately for each constellation while our nML 
detector is able to detect arbitrary constellations without any 
modification. The mixed use of one-bit and high-resolution 
ADCs in massive MIMO was analyzed in |[30l, where it was 


found that using one-bit ADCs with a few high-resolution 
ADCs can achieve similar performance with the unquantized 
case. Joint channel and data estimation using low-resolution 
ADCs was proposed in ED- While the performance was 
comparable to the case with perfect CSI, the computational 
complexity of the joint technique OTI may still be too high 
to be affordable in a commercial system. 

Since our detector exploits CSI as part of the detection, we 
also propose an ML channel estimation technique using one- 
bit ADCs that is in line with our proposed nML detectors. 
There has been related work on channel estimation with low- 
resolution ADCs ll^ and with one-bit ADCs E3-E21- In 
ll34l . the expectation-maximization (LM) algorithm was 
exploited to estimate channels. The problem solved in ES, 
ED, however, is not convex, and the LM algorithm may 
converge to a local optimal in the high SNR regime. The 
work 051 . used generalized approximate message passing 
(GAMP) algorithms for channel estimation with one-bit ADCs 
where the techniques heavily relied on the sparsity of the target 
vector. While the work in El, El was not able to estimate 
the norm of the channel due to the simple zero-threshold 
setting for one-bit quantization, ESl adopted asymmetric one- 
bit quantizers to deliver the the norm information of the target 
vector. Our ML channel estimator, which is an extension of 
El, is shown to estimate not only the channel direction but 
also the channel norm and does not make an assumption about 
sparsity in the channel. 

Our contributions are summarized as follows. 

• We propose an nML detector for uplink multiuser massive 
MIMO systems. The proposed nML detector is based on 
the ML detector developed for distributed reception in 
ED- We show that the two problems, i.e., distributed 
reception and multiuser detection, are essentially the 
same problem, which makes it possible to exploit the 
detectors from ESI for our multiuser detection problem. 
The complexity of the ML detector in ESl . however, 
grows exponentially with the number of uplink users, 
which prevents its use in massive MIMO with many 
users. We reformulate the ML detector in ESl to derive 
an ML estimator and convert the ML estimation problem 
into a convex problem. Therefore, we can rely on efficient 
convex optimization techniques to obtain an ML estimate 
and perform symbol-by-symbol detection based on the 
ML estimate. 

• We implement a two-stage nML detector to further im¬ 
prove the detection performance. The two-stage nML 
detector exploits the structure of the original ML detector 
with the reduced candidate set constructed by the pro¬ 
posed (one-stage) nML detector. Numerical results show 
that the two-stage nML detector improves the detection 
performance significantly in the high SNR regime and 
achieves similar performance with the original ML de¬ 
tector. 

• We derive the exact probability that two different transmit 
vectors (after they have passed through the channel and 
noise is added) result in the same quantized received 
signal with one-bit ADCs in a certain condition. The 
result shows the relationship between the probability and 
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the numbers of antennas and users, where the probability 
of having the same quantized outputs goes to zero as the 
number of antennas goes to infinity. 

• We propose an ML channel estimator that has the same 
structure with the proposed nML detector. The proposed 
estimator can estimate the direction and norm of the 
channel more accurately than other channel estimators 
using one-bit ADCs. Because of the similar structure, it 
is possible to implement both the nML detector and the 
ML channel estimator using the same algorithm. 

The paper is organized as follows. We describe our system 
model using one-bit ADCs in Section |II] In Section |III] we 
briefly discuss the detectors which were originally developed 
for distributed reception in Then, we propose our nML 
detectors and present the asymptotic analyses in Section |IV] 
We propose an ML channel estimator in Section |V] In Sec¬ 
tion m we evaluate the proposed techniques by simulations, 
and the conclusion follows in Section fVlIl 

Notation: Lower and upper boldface letters represent col¬ 
umn vectors and matrices, respectively. ||a|| is used to denote 
the £ 2 -norm of a vector a, and A^, A*, A^ denote the 
transpose, Hermitian transpose, and pseudo inverse of the 
matrix A, respectively. Re(b) and Im(b) represent the real 
and complex part of a complex vector b, respectively. Om is 
used for the m x 1 all zero vector, and Im denotes the mxm 
identity matrix. and represent the set of all m x n 

complex and real matrices, respectively. 

11. System Model 

We explain our system model and several assumptions that 
are relevant to our detector design. We also define expressions 
for one-bit ADCs in this section. 

A. Massive MIMO Received Signal Model 

We consider a uplink multiuser cellular system with 
cells. Each cell consists of a base station with Aj. received 
antennas and K users equipped with a single transmit antenna. 
All KNc users transmit independent data symbols simulta¬ 
neously to their serving base stations. Assuming all users 
transmit data with power P, the received signal at the i-th 
base station = [yi^i y 2 ,i • • • yNr.i] is 

K Nc K 

y^ = v^E v^EE Hz (1) 

k—1 m—1 k—1 

m^i 


such that 

llx,||" = A (2) 

for all i where = \xii xa ■ ■ ■ XiK^■ If 5 is a phase 
shift keying (PSK) constellation satisfying |smP = 1 for all 
m, then the norm constraint is trivially satisfied. For a QAM 
constellation, the law of large numbers (LLN) gives 

K 

Y,\x,k?-K (3) 

fc=i 

as K becomes large if all users select their data symbols 
independently and with equal probability within S (with 
proper normalization). Although the norm constraint in da 
plays an important role in implementing the proposed nML 
detector, numerical results in Section |VT] shows that the nML 
detector works well even with moderate numbers of users. 

Assumption 2: We assume that each base station has perfect 
local CSI with which to implement its detector. After imple¬ 
menting the detectors, we relax this assumption in Section IV] 
where we consider channel estimation techniques for one- 
bit ADCs. We neglect pilot contamination during the channel 
estimation procedure because the channel is already severely 
distorted due to one-bit ADCs. Note that pilot contamination, 
which contaminates the channel estimate, remains the only 
channel impairment with full-resolution ADCs and extremely 
large number of antennas lEH, iQl. 

In addition to these assumptions, we first focus on the single 
cell scenario because the detectors considered in this paper 
do not exploit any kind of inter-cell cooperation. In Section 
IIV-DI we explain how the proposed detector can be adapted 
to a multicell setting. 

For the single cell scenario, the received signal y in ([T]) 
becomes 

K 

y = V^E hfcTfc-f n = V^Hx-f n, (4) 

fc=i 

and the SNR is 

P 

P=—- (5) 

We assume the base station knows p perfectly. 

B. Received Signal Representation with One-Bit ADCs 

We focus on a massive MIMO system that uses one-bit 
ADCs for the real and imaginary parts of the each element of 
y. The conceptual figure of our system is depicted in Fig. [T] 
The output of the n-th receive antenna after the one-bit 
ADCs is given as 


where h.i^rnk & is the channel vector between the i-th 

base station and the fc-th user associated with the m-th base 
station, Xmk is the data symbol, which satisfies E [x^k] = 0 
and E [|a:mfcP] = 1, transmitted from the fc-th user supported 
by the m-th base station, and rii ~ CA/'( 07 Vr,cr^IjvJ is the 
complex additive white Gaussian noise (AWGN) at the i-th 
base station. 

Before developing the detectors, we make the following 
assumptions. 

Assumption 1: The data symbols Xik are from an M-ary 
constellation S = {si,--- ,Sm} and have been normalized 


ijn = sgn(Re(t/„)) -f jsgn(Im(?/„)) 
where sgn(-) is the sign function which is defined as 

1 if X > 0 


sgn(x) = 


— 1 if X < 0 


Therefore, we have 

Vu G {1 +J,-1 +J,-1 - j, 1 - j} 
for 1 < n < Ap. The collection of ijn is given as 

y=[yi 2/2 ••• 2 /tvE- 


( 6 ) 

(7) 

( 8 ) 

(9) 
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K users 


Fig. 1: MU-MIMO with K users and receive antennas. 
Each received signal is processed with two one-bit ADCs. 


III. Possible Detectors Using One-Bit ADCs 

In this section, we reformulate the detectors for distributed 
reception proposed in IITSl into the uplink multiuser massive 
MIMO setting. We also discuss their characteristics and limi¬ 
tations. This discussion is useful for developing our detectors 
in Section HVl 


A. ML Detector Reformulation 

Let G be the n-th row of the channel matrix H, 

i.e., 


H = gl g2 


gN,. 


( 10 ) 


Note that g„ is the channel between the K users and the n-th 
receive antenna. Assuming the real and imaginary components 
of the Gaussian noise are IID, it is useful when approaching 
one-bit ADC problems to rewrite the signal model in the real 
vector form instead of complex form as 


— 




XR = 


Re(g„) Im(g„) 
-Im(g„) Re(g„) 


Re(n„) 

Im(n„) 

Re(x 

Im(x 


1 T 

gR.n.l’ 


.gR.".2. 


I, 2 x 2 ir 



I^R,n,l 


_'RR,n,2_ 


i,2R:x1 


where 


and 


SR,n,l — 


R.e(g„) 


Im(g„) 

_-Im(gn) 

) SR,n,2 — 

Re(g„)_ 


^R,n,2 


v|o,;l 


( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 


for all n and i. The received signal at the n-th receive antenna 
can be also rewritten as 


yR,n — 


R-e(y„) 


yR.n.l 

_Im( 2 /„) 


yR,n.2 


= VpGr,„xr -f nR,„, (16) 


and the vectorized version of the quantized jjk in the real 
domain is given as 


yR,n — 


sgn(Re(y„)) 


yR,n,l 

_sgn(Im(y„))_ 


yR,n,2_ 


(17) 


Based on yR^„, the base station generates the sign-refined 
channel matrix for the n-th receive antenna as 


SR,n,l 
.gR.n,2. 

where gR,„,i is defined as 



gR,n,z — yR,n,zgR,n ,2 • 


(18) 


(19) 


Note that yR,„,i = ±1 depending on the sign of yB.,n,i- Dehne 
5 r to be 


5r 


Re(si) 

Im(si) 


Re(sM) 

Im(sM) 


( 20 ) 


where M is the size of the data symbol constellation S. 

With these definitions and using similar logic from 1^ . the 
ML detector can be defined as 


xr.ml = argmax 


n n $ (v^gR,„,iXR) 


( 21 ) 


where $(f) = dr and is the AT-ary Carte¬ 

sian product set of 5r, which is ordered with the real parts of 
the constellations first and the imaginary parts later. The \/2 
term in dSlT i comes from the distribution of nn^n,i given in 


B. ZF-Type Detector Reformulation 

By brute-force search, the complexity of the ML detector 
in (ISTT l is which grows exponentially with the number 
of users K. To support a large number of K, the ZL-type 
detector was proposed in ll3^ . The base station first obtains 
the ZL estimate as 

xzF = Hty. (22) 

Because the norm square of xzp may not equal to K, the base 
station normalizes Xzf as 

XZF = (23) 

IIxzfII 

and performs symbol-by-symbol detection using xzf as 

TzF.fe = argmin \xzF,k - xf (24) 

where xzF,k is the fc-th element of normalized xzf- The 
normalization is not an issue for PSK constellations; however, 
it is crucial for QAM constellations. 

It was shown numerically in ll^ that ZL-type detectors 
saturate at a higher error rate floor than ML detectors as SNR 
increases; the error rate floors in both cases are inevitable due 
to the one-bit ADCs. Therefore, we propose new detectors that 
outperform the ZL-type detector and require less complexity 
than the ML detector. 
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IV. Near ML Detector Implementation 

We derive our nML detector by converting the original ML 
detection problem as a convex optimization problem. This can 
be done by relaxing constraints on the transmitted vector. We 
also analyze the nML detector in the asymptotic regimes and 
derive the exact probability that two different transmit vectors 
result in the same quantized signal (after they have passed 
through the channel and noise is added) in a certain condition. 
The relationship between the probability and the number of 
receive antennas show that the probability goes to zero as the 
number of antennas goes to infinity, which shows the benefit 
of massive MIMO for one-bit ADCs. Then we implement 
the two-stage nML detector to further improve the detection 
performance. Finally, we extend the proposed detectors to a 
multicell scenario. 


A. Convex Optimization Formulation of nML Detector 

Because of the norm constraint 0, we define the ML 
estimator by relaxing the constraint xr £ in (|2T]) as 

2 Nr 

^R.iiL = argmax 11 ® (v^gR.n.zXR) (25) 

Uk\\^=K 

2 Nr 

= argmax V V log $ . (26) 

xrEM X n=l 

i|xr||"=a 

It was shown in Il38l that 

^R,ML (27) 

in probability as iVr —^ cxj for arbitrary p > 0. Therefore, if 
(|2^ can be solved, then the detector is guaranteed to achieve 
good performance when N,- is large. In fact, all of the properly 
designed detectors including linear detectors become optimal 
as TVr —>■ oo m, m, and (l27l i proves that the ML detector 
based on (l26l l is indeed a proper detector. 

The maximization in (l26l l is still not easy to solve in general. 
The function $(•) is log-concave but the optimization problem 
in (|2^ is not convex due to the norm constraint ||xr|P = K. 
To sidestep this challenge, we relax the norm constraint as 
||xr|P < K and reformulate the problem as 

2 Nr 

Xr,ml = argmax X! X! ^ (28) 

xrGR"^><i i=i „=i ^ 2 

I|xr|F<k 


which is a convex optimization problem that can be efficiently 
solved nm. Similar to Il29l , 14^ . we provide a simple yet 
effective iterative approach to solve (l28l l in Algorithm [T] 

In Step 1, we initialize to 


^(0) 


= '/k 


Gjl 


R^2JV, 


GSl 


rJ-2A, 


(29) 


to resemble the maximum ratio estimate. The quantization 
using one-bit ADCs has been reflected in Gr as shown in 
(US, therefore, the maximum ratio estimate is based on the 
all-one vector. Steps 6 and 7 are based on projected gradient 


Algorithm 1 Iterative algorithm to solve 


,(o) 


Initialization 

1; Set the initial point x^ 

2: Set 

Gr = [gR,l,l ••• gR.LAr gR.2,1 ■■■ gR,2,AR] 

3; Set the step size k and the termination threshold e 

Iterative update 


4: 

5; 


6 : 

7 

8 : 

9; 


While 


dk) 


.(k-1) 


(fc) _ Nk-l) 

T3 - -^T3 


> e 


■nGl 


V/ 


‘r 


where the i-th element of V/ (z) is 

- 1 


If 


'(fe) 


■R 
,{k) 

A 

end If 
end While 


> K 

= Vk 




Xr 


method ll43l to ensure the norm constraint. Note that a similar 
iterative algorithm was proposed in 1291 . While our approach 
can support arbitrary constellations without any modification 
on Step 6 and 7, the constraint on the algorithm in l2^ is 
separately adapted to each constellation, i.e., for M-QAM 
constellation, the constraint on each element of xj ^' is set to 
1)A, (V^M — 1)A] in l29l where A is the average 
power normalization factor. 

( 2 ) 

After obtaining the estimate xj^ the base station needs to 
perform normalization followed by symbol-by-symbol detec¬ 
tion similar to the ZF-type detector in (l24l l. If we let xr^ml./c 
be the fc-th element of 


Xr,mL 


= Vk 


J2) 

%,ML 


^(2) 

xr,ML 


(30) 


the nML symbol-by-symbol detection is 


^nML.fc = argmin |(a;R,ML.fc + Jxr,ml,K+ k) - x\ (31) 

x^S 

considering the fact that xr^ml is a 2K x 1 real vector. 

Remark 1: As also discussed in 1291 . the complexity of 
Algorithm 1 is dominated by the number of iterations because 
there are only two matrix-vector multiplications (one in the 
initialization in 12^ and the other in the iteration loop in 
Algorithm 1) and numerical calculations to get V/(z)j in 
the iterative update. The numerical calculations for updat¬ 
ing Vf{z)- in each iteration prevent performing the direct 
comparison with the ZF-type detector but it is reasonable to 
assume that the complexity of one iteration is less than that of 
the ZF-type detector that requires matrix inversion. As shown 
in Section IVII the number of iterations for Algorithm 1 to 
converge can be moderate, e.g., less than 20. 

Remark 2: Note that either nML detector based on (|2^ or 
(l28l is suboptimal compared to the original ML detector (l2Tll . 
The proposed nML detector is based on the ML estimation 
(l28l l over the 2iT-dimensional real space with the norm 
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constraint. It can be the case that 


2 Nr 

^ ^ log$ 

i—1 n—1 

2 Nr 

> (\/^gR,„ ,zXr) 

2=1 n—1 
2 Nr 

> X! H (\/^gR,n,*XR) 

2=1 n —1 


(32) 

(33) 

(34) 


with 


.^(2) 

2 

-(2) 


> 

^R,ML 


(35) 


where xr is the true transmitted vector and xr S \ {xr}. 
In words, it can happen that the original ML detector makes 
the correct decision as in (l33l l and (l34l l while the estimated 
vector from (l28T l is closer to xr that is not the true 

transmitted vector. The suboptimality, however, will not come 
into play when is large as shown in |[38l . 


B. Analyses in Asymptotic Regimes 

The estimates in (l26l l and in (l28T l may not be 

the same in general because is selected from a circle 

||xr|P = K while xj^M^ is selected from a ball ||xr|P < K. 
Note that (|26] | can have multiple local optimal solutions due 
to non-convexity of the constraint. If we let be a set of 
all possible solutions of (l26l l. the following lemma shows the 
relation between and xj^M^ in the high SNR regime. 

Lemma 1. When p ^ oo, we have 

(36) 

Proof: Please see Appendix lAl ■ 

Lemma [U shows that the nML detector based on (l28l l will 
perform the same as in (l26T l with the relaxed norm constraint 
in the high SNR regime. 

It is necessary that different transmit vectors result in 
different quantized received signals to avoid possible detection 
errors because gR,n,i in (fT9l l is based on the sign refinement 
from the quantized received signal y’s.,n,i- Although it is 
difficult in general, we can explicitly derive the probability 
of which two different transmit vectors result in the same 
quantized received signal for a special case. 

Special case; Consider two transmit vectors xi = 
[a;i X 2 ■ ■ ■ xk\ and X 2 = [—xi X 2 ■ ■ ■ Xk\ 

where Xk, which is selected from a standard M-ary constel¬ 
lation S with equal probability, is the transmit symbol of the 
A:-th user. The two received signals are yi = v^Hxi -|- n 
and y 2 = V^Hx 2 -f n. Define yj^ and y 2 as the quantized 
outputs of yi and y 2 using one-bit ADCs. Assume each entry 
of H follows IID Rayleigh fading and n ^ CAf(OAr,., 

Proposition 1. For the special case, 

Pr (yi = y 2 ) = ^ arctan . (37) 


The following corollary is a direct consequence of Propo¬ 
sition [T] 

Corollary 1. For the special case, Pr (yi = y 2 ) —>■ 0 as 
N-[ —7> oo. 

Corollary [T] shows that, for arbitrary SNR values, the two 
transmit vectors from the special case give different quantized 
signals as the number of receive antennas goes to infinity, 
which shows the benefit of massive MIMO with one-bit ADCs. 
Although it is hard to generalize Corollary [T] for arbitrary pair 
of transmit vectors, we expect a similar result would hold in 
general. 


C. Two-Stage nML Detector 


To improve the performance of nML, we also propose a 
two-stage nML detector. The two-stage nML detector reduces 
the number of candidate transmit vectors, based on the output 
of the one-stage nML detector. Using the estimate xr_ml in 
(l30l l and the detected symbol XnMh,k in OTl l. we define the 
candidate set of the /c-th element 

|(^R,ML,fc + j^R.ML.if-rfc) — x\ 

|(a:R,ML,/c + jXR,ML,K-Vk) — ^nML.fe 

(38) 

where c is a constant that controls the size of Xk- With Xk 
for k = 1,..., K, we can define the reduced candidate set of 
the transmit vectors as 





If we let 


A’r = 


Re (x) Im (x) 


x€X 


(39) 

(40) 


then the second stage of the two-stage nML detector can be 
defined as 

2 Nr 

XR,two-nML = argmax nn 4* (v^gR.n.iXRj . (41) 

The set X (or Ar) allows testing possible transmit vectors 
that are close to the estimate xr uML- If we let the constant 
c — 00 , the two-stage nML detector becomes the original 
ML detector. With a proper value of c, the numerical results 
in Section|VT]show that it is possible to effectively improve the 
detection performance especially in the high SNR regime with 
marginal additional computational complexity. The concept 
of the two-stage nML detector is similar to that of sphere 
decoding that exploits reduced search space. 


D. Extension to Multicell Setting 

For the multicell scenario, the inter-cell interference should 
be taken into account for the proposed detectors. It is reason¬ 
able to assume that the base station can accurately estimate 
the long-term statistic of the inter-cell interference. Without 
having any instantaneous CSI from out-of-cell users, however, 
the base station will consider the inter-cell interference as 
additional AWGN, i.e., 

Nc K 

v^EE (42) 

m=l k^l 
m^i 


Proof: Please see Appendix iBl 
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where rji captures the long-term inter-cell interference statistic, 
e.g., pathloss and shadowing. Then the effective signal-to- 
interference-noise ratio (SINR) at the i-th base station is 


Pi,MC 


P 

Pr]f + cr^ ’ 


(43) 


and the proposed detectors can be adapted to the multicell 
scenario by substituting p in (l28l l and (l4Tl l with pi,MC- 


V. Channel Estimation with One-Bit ADCs 

For a coherent detector, the CSI is normally obtained 
through an estimate of the channel. In this section, we develop 
channel estimation techniques for one-bit ADCs at the base 
station. We focus on estimating g„, i.e., the channel between 
the receive antenna n and K users, instead of h^. 

We consider a block fading channel to develop channel 
estimation techniques. We assume the channel is static for 
L channel uses in a given fading block and changes indepen¬ 
dently from block-to-block. The received signal at the n-th 
antenna for the f-th channel use in the m-th fading block is 
given as 


4“ (44) 

Let the first T < L channel uses be devoted for a training 
phase and the remaining L — T channel uses be dedicated to 
a data communication phase. Put the first T received signals 
during the training phase into a vector form as 


yn,m, train — 

,trainSn,m “t" '^n,m,train 

(45) 

where 



yn,m, train ~ [yn,m[0] 

yn,m[T-l]]* 

(46) 

^m,train — [^m [0] * ’ ’ 

Xm[T-l]] GC^^^, 

(47) 

W7T,^y7j^train — [^n,m[0] 

- Wn,m[T - 1]]* G 

(48) 

In the training phase, Xm^train is known to the base station 


but gn,m must be estimated. While arbitrary training matrices 
are possible, for simulation purpose in Section IVII we focus 
on unitary training where train satisfies ii 

X™,trainX:„,train = TIk if K < T. (49) 


The normalization term T ensures the average transmit power 
equals to P in each channel use. 

Similar to the previous sections, we reformulate all expres¬ 
sions into the real domain as 


yR,n,m,train — "x/P^Rim.trainSR.njm 4“ train (50) 


where 


yR ,n,m,train 


X 


R,m,train 


SR,n,m 


Wr 

,n,m,train 


(yn,m,train 
Im (yn,??i, train 

R^(X^ train) 
Im(Xyj.j train) 


p2Txl 


Im(Xyj.j train) 
Rl^(X^ train) 


(Sn,m) 

Im 


e 


R,e train) 

Im train) 


G R2 Tx 1 


(51) 

(52) 

(53) 

(54) 


It is important to point out that (fSOl l has the same form as 
(fTbl l while the roles of the channel and the transmitted signal 
are reversed. Therefore, using the same techniques that we 
exploited for the detectors, we can develop channel estimators 
based on the one-bit ADC outputs and Xr ^, train- 

We define the f-th column of Xr,„, train as train,i and 

the i-th output of the one-bit ADC as 


yR,n,m,train,i — Sgll(yjin,m,train,i) (55) 

where ytl^n, 7 n,train,i is the i-th element of yR,„,m, train- Note 
that there are 2T one-bit ADC outputs in total for the n- 
th receive antenna, i.e., T outputs for each of the real and 
imaginary parts of the received signal. Based on yR,k,m,train,i, 
the base station performs the sign-refinement as 


^R,m,train,i — yR,n,rn,train, 0 ^R,m,train,i-/ 


(56) 


and the ML channel estimator is given as 

2T 


gR.n,m.ML = urgmax 
gRGR2^ 


n ^ (v^^R ,m,train,i Sr) (57) 


2T 


= argmax V log (t> ( ^ ,gR 


(58) 


Because $(-) is a log-concave function, and there is no 
constraint on gR, it is possible to solve (fSST l using standard 
convex optimization methods (or suitably modified Algorithm 
1). Unfortunately, this channel estimator tends to overestimate 
the norm 


IISr, n,m, ML II ^ IISR, n,m II (59) 

due to the fact that log$(-) is an increasing function where 
SR,n,m is the true channel vector. 

To overcome this problem, we impose a norm constraint on 
gR and convert (fSSl) to 


2T 

gR,n,m,ML = argmax V log (t> ( V^Xr 
ErGR^'^xi,^^ V V 

l|gRll"<-R 


train,! gR 


(60) 


using the fact that E ||gR,„,m|| = K for most of channel 

models. The norm constraint can be further optimized if the 
base station knows the long-term statistic of the channel norm. 

For comparison purposes, we also define a simple ZF-type 
channel estimator as 


gR,n,m,ZF 


= y/K 


^R,m,train 


)^yR ,n,m,train 


(^R,m.train) yR,n,m,train 


(61) 


which is forced to satisfy ||gR,n,m,ZF|P = K. We numerically 
compare our ML channel estimator and the ZF-type channel 
estimator in Section |VT] 

















SNR (dB) 


Fig. 2: SER vs. SNR (in dB) of four detectors with = 32, 
K — 4, and M — 4 (QPSK). The original ML, the one- 
stage nML, the two-stage nML, and the ZF-type detectors are 
compared. 



Fig. 3: SER vs. SNR in dB scale with M = 8 (8PSK), Ff = 8, 
and different values of N-^. The GAMP detector is from ll28l . 


VI. Simulation Results 

We perform Monte-Carlo simulations to evaluate the pro¬ 
posed techniques. The one-stage nML detector is based on 
Algorithm 1 with a termination threshold e = 10“^ and step 
size K = 0.01. For the two-stage nML detector, we let the 
constant c = 1.3 in (l38] l. We first consider the single cell 
scenario to compare detector performance. Then we take the 
multicell scenario into account. 

A. Single Cell Scenario 

We assume IID Rayleigh fading channels, i.e., all elements 
of H are distributed as CA/^(0,1), although the distribution of 
the channel was not explicitly incorporated into the proposed 
detector designs. For the time being, we assume the base 
station has perfect CSI and evaluate the detectors. Later, we 



Fig. 4: SER vs. SNR in dB scale with M = 16 (16QAM), 
K — 8, and different values of Wj.. The GAMP detector is 
from ll28l . 


evaluate the detectors with imperfect CSI. We use the average 
symbol error rate (SER) which is defined as 

1 ^ 

SER = {Xn ^ Xn I X Sent, H, n, p, K, N^, 5)] 

n—1 

(62) 

for the performance metric where the expectation is taken over 
X, H, and n. 

We first compare the four detectors: 1) the original ML 
detector (1211 1 that is based on exhaustive search over all 
possible transmitted vectors, 2) the one-stage nML detector 
(l28l l that is based on convex optimization, 3) the two-stage 
nML detector explained in Section IIV-CI and 4) the ZF-type 
detector. Due to the computational complexity of the original 
ML detector, we set AT = 4, M = 4 (QPSK) for all users, and 
Wr = 32 which may not be considered as massive MIMO. 
We plot SERs of the four detectors in Fig. |2] The figure 
shows that the one-stage nML detector is suboptimal compared 
to the original ML detector as discussed in Remark 2 of 
Section IIV-AI while the two-stage nML detector gives almost 
the same performance with the original ML detector. Further, 
both the nML detectors outperform the ZF-type detector. 

In Fig. [3 we plot the SERs of the one- and two-stage nML 
detectors, the ZF-type detector, and the detector based on the 
GAMP algorithm from 1^ according to SNR with K = 8, 
M = 8 (8PSK), and different values of N,-. We have modihed 
the GAMP algorithm to have adaptive step sizes as in m 
for better performance. We do not consider the original ML 
detector in this scenario because of its excessive complexity, 
i.e., the detector needs to compare 8® = 2^^ possible transmit 
vectors. The figure shows that the two-stage nML detector 
and GAMP detector are comparable and outperform the one- 
stage nML and ZF-type detectors with the same number of 
N,-. The ZF-type detector suffers from an error rate floor while 
other detectors do not have such floor until 10“® SER. With 
small SERs, it is possible to use weaker channel coding (with 
higher code rate or shorter block length) to improve the system 
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Fig. 5: MSE and normalized MSE of different channel estimators with K = 8 and different values of T and SNR . Expectation- 
maximization (EM) estimation is from 0^ . and “EM, fixed norm” is the same as EM except the norm square fixed to K. 


TABLE I; Average number of iterations for Algorithm 1 and 
size of X in (|39] | with K = 8. 



Algorithm 1 

Size of X 

M = S, Nr = 128 

18.3721 

2.6475 

M = S,Nr = 192 

12.1172 

1.9522 

M = 16, Nr = 128 

18.8378 

4.7196 

M = 16, Nr = 192 

13.0174 

4.6471 


throughput, which has been shown using binary symmetric 
channels in ll46ll . 

We adopt the same system setup with Eig. [3 except a M = 
16 (16QAM) constellation for data symbol in Eig.|4] The two- 
stage nML detector and GAMP detector are still comparable; 
however, the two detectors also suffer from an error rate floor. 
With = 196, the error rate floors are mitigated for all 
detectors, which shows the benefit of massive MIMO for using 
one-bit ADCs. 

Eigs. [3 and |4] both show that the proposed detectors give 
much better SER performance compared to the ZE-type detec¬ 
tor with the same number of receive antennas, which means 
that the proposed detectors are able to support more users than 
the ZE-type detector with the same system setup. Therefore, 
if the base station has sufficient computational power, it is 
always beneficial to use the proposed detectors than the ZE- 
type detector. Note that the computational power at the base 
station is related to the digital baseband processing in Eig. 
[T] and different from having one-bit ADCs. Therefore, the 
benefit of using one-bit ADCs, e.g., power consumption and 
cost, still holds for the proposed detectors although higher 
computational power is required for the digital baseband 
processing at the base station. 

To evaluate the complexity of the proposed detectors, we 
compare the average number of iterations for Algorithm 1 to 
converge for the one-stage nML detector and the average size 
of X in for the two-stage nML detector, with K = 8 
and different values of M and in Table |I1 It shows that 


Algorithm 1 requires less than 20 iterations to converge in 
average and the additional comparison, i.e., the size of X, 
performed in the two-stage nML detector is marginal. It is 
interesting that the iteration number of Algorithm 1 decreases 
with the number of antennas, which shows that the proposed 
detectors based on Algorithm 1 will become more efficient 
with large Wj.. 

Now, we evaluate the ML and ZE-type channel estimators 
discussed in Section|V] We focus on estimating gn,m, i-S-, the 
channel between the n-th receive antenna and K users. We 
define the mean squared error (MSE) of a channel estimator 
X as 


1 


MSEx = -E 


|Sn,m Sn,m,x\ 

and the normalized MSE (NMSE) as 


NMSEx = -^E 
K 


gn 


g 


n,m^x 


(63) 


(64) 


which are used as performance metrics. The expectations are 
taken over gn,m- In Eig.|5] we compare the proposed ML chan¬ 
nel estimator to the ZE-type estimator and the expectation- 
maximization (EM) method from ll34l with different training 
lengths T and SNR values with K = 8. Regarding MSE, the 
proposed ML estimator outperforms other estimators for both 
p = 0 and 20dB cases. The EM method performs well only 
when T is large with p = OdB and fails to estimate the channel 
norm when p = 20dB as also shown in ll34l . The EM method 
with the norm fixed to K gives better performance than the 
ZE-type estimator when p — OdB while the two estimators 
become exactly the same with p = 20dB. Regarding NMSE, 
the proposed ML estimator and the EM method (and the one 
with the fixed norm as well) are comparable when p = OdB 
while the ML estimator outperforms the EM method with 
p = 20dB. 

To verify the effect of channel estimation, we plot the 
SERs of the two-stage nML detector and GAMP detector with 
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SNR (dB) 


Fig. 6; SER vs. SNR in dB scale with M = 16 (16QAM), 
K — 8, and — 192. The two-stage nML detector and 
the GAMP detector from ll28l are compared with different 
assumptions on channel information. 

different assumptions on CSI; perfect CSl, perfect GDI, and 
GDI with error that is defined as 

Sn,m/||Sn,m|| T V^NMSE© 7 ^^yj.j/|[© 7 j^m|| /gTc\ 

-^- (65) 

Sn,m/\\Sn,m\\ + VNMSE©„,m/||©„,m|| 

where the element of ©„,m is distributed as CAf{0, 1). We set 
K = 8, M = 1Q, N, = 192, and NMSE = IQ-^. It is shown 
in Fig. |6] that the two-stage nML detector outperforms the 
GAMP detector using GDIs when SNR is low, which shows 
that the GAMP detector requires more accurate channel norm 
information than the proposed detector for low SNR values. 

To compare the detectors in a practical setting, we combine 
the detectors with a low-density-parity-check (LDPG) code. 
We assume the base station has perfect GSI for this study. We 
adopt a rate 1/2 LDPG code with the block length of 672 bits 
from the IEEE 802.1 lad standard BtI . After hard detection 
by the detectors, the estimated symbols (or bits) are decoded 
using the bit-flipping decoding algorithm ll48l . The coded bit 
error rates (BERs) of the one-stage nML detector and ZF- 
type detector according to SNR with K = A, N = 64, and 
M = 8 are shown in Fig. |7] The figure clearly shows that the 
nML detector outperforms the ZF-type detector even for this 
practical setting. Further improvements could be expected if 
further work is put into deriving an appropriate soft decision 
decoding metric, which requires the probability distribution of 
ML ll28l l. It is also possible to exploit the approximated 
soft metric as in Il49l . but we leave this for future work. 

B. Multicell Scenario 

For the last numerical study, we consider the multicell 
scenario where the detailed simulation parameters are listed in 
Table Hn We consider two different user dropping scenarios. 1) 
All users except a typical user in the center cell are randomly 
dropped within corresponding cells. The typical user is located 
with the distance d (and a random angle per drop) from the 



SNR (dB) 


Fig. 7; Goded BER vs. SNR in dB scale with M = 8 (8PSK), 
K — A, and = 64. The proposed nML and ZE-type 
detectors are compared using the rate 1/2 LDPG code adopted 
in the IEEE 802.1 lad standard. 


base station in the center cell. 2) All users except the users 
in the center cell are randomly dropped within corresponding 
cells. The users in the center cell are randomly dropped within 
the range (d—20m, d-|-20m) from the center cell base station. 
The second scenario can be considered as coordinated uplink 
transmission with user scheduling that selects users with sim¬ 
ilar received signal power while the first scenario corresponds 
to uncoordinated uplink transmission. We consider coded BER 
of the typical user for the hrst scenario while BERs of all K 
users in the center cell are averaged for the second scenario. 
We neglect shadowing effect and rji in (l42l i is based on the 
distances between the center cell base station and out-of-cell 
users. 

The BER results for these two scenarios according the 
distance d are plotted in Eig. [8] We can see that the one-stage 
nML detector outperforms the ZF-type detector for both user 
dropping scenarios. As d increases, the BER performance of 
both detectors becomes worse because of the reduced received 
signal power. Note that the BER performance of the second 
scenario is much better than that of the hrst scenario when 
d is large. Eor large d, the received signal power of the 
typical user in the hrst scenario is overwhelmed by other 
users’ received signals (the near-far effect), resulting in poor 
BER performance. If all users experience similar SINR as 
in the second scenario, however, the BER performance is 
quite good even with large d. This shows that the nML 
detector will perform well with proper user scheduling or 
uplink power control, which are already common in current 
cellular systems. Note that the ZF-type detector suffers from 
the notable error rate hoor for the second scenario when d 
is small (that corresponds to the high SNR regime for the 
single cell scenario) while the proposed nML detector does not 
have such floor until 10“® BER. In the first scenario, there is 
no error rate floor even for the ZE-type detector because the 
received signal of the typical user overwhelms other users’ 
received signals. 
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TABLE II: Multicell Simulation Parameters 


Parameter Assumption 


Cell layout 
Cell radius 

# of RX antennas per BS (A^r) 
# of TX antennas per user 
# of users per cell (K) 
Min. dist. btw. BS and user 
User transmit power 
Path loss per km 
System bandwidth 
Noise spectral density 
Noise figure 
Constellation 
Channel coding 


57 hexagonal cells 
500m 
64 
1 
4 

100m 

23dBm 

131.1+42.81ogio(dist.) dB 
5MHz 

-174dBm/Hz 

5dB 

8PSK 

Rate 1/2 LDPC from 802.1 lad 



d(m) 


Fig. 8; Coded BER vs. d (in meters) for the multicell setting 
with parameters in Table HH Coordinated and uncoordinated 
user dropping scenarios are compared for the nML and ZF- 
type detectors. 


same hardware can be used for both tasks, which make the 
proposed techniques more attractive for uplink massive MIMO 
systems using one-bit ADCs. 

There are several directions for future research. The pro¬ 
posed detectors can be extended to frequency selective channel 
considering their possible use in mmWave systems, for exam¬ 
ple by extending the results in ll25l . ifSOl . To make the proposed 
ML channel estimator more practical, the training overhead 
should be reduced. As in ED, El, it may be possible 
for the proposed ML channel estimator to exploit channel 
statistics, e.g., the temporal and spatial correlation, to reduce 
the training overhead. It would be also interesting to consider 
maximum a posteriori (MAP) detector and channel estimator 
using one-bit ADCs and compare them with the proposed 
techniques where MAP detector and channel estimator with 
low-resolution ADCs have been studied in EOl . 


Appendix A 
Proof of LemmaEI 


Recall the estimator in (l281 l. When p —?> oo. 


oo 


( 66 ) 


unless ^ jXR = 0. Therefore, the estimator ( |28] ) hnds 

(2) ’ ’ 

^R ML satisfies 


gR.n,zxgl^L > 0 (67) 

for all 1 < n < W]. and 1 < i < 2 in the high SNR regime 
because $(f) is an increasing function of t but upper bounded 
by 1. Then, 


VII. Conclusion 

We proposed a unified framework of detection and channel 
estimation techniques for uplink multiuser massive MIMO 
systems using a pair of one-bit ADCs at each antenna. 
The proposed techniques are based on off-the-shelf convex 
optimization methods, which makes it easy to implement in 
practice. We proposed two nML detectors, i.e., one-stage and 
two-stage nML detectors, which give better performance than 
the ZF-type detector for all range of SNR regimes and number 
of antennas. The two-stage nML detector reduces the number 
of candidate transmit vectors using the output of the one-stage 
nML detector and exploits the ML detector structure to further 
improve detection performance of the one-stage nML detector. 
Numerical studies showed that the proposed detectors are able 
to perform well even with not-so-large number of antennas, 
robust to inaccurate channel estimation, and outperform the 
ZF-type detector for practical channel coded and multicell 
settings. We also proposed a ML channel estimator that can 
effectively estimate not only the direction but also the norm of 
the channel even with one-bit ADCs. Because of the unihed 
structure of the proposed detectors and channel estimator. 


2 Nr 

^ ^ log $ ( V^gR,„.ixg^ML) (68) 

i—l n—1 

2 Nr 

> X] S ^ {V^SR,n,ia^RMh) (69) 

i—l n—1 

with arbitrary 0 < a < 1. Therefore, the norm square of 
( 2 ) 

^R ML always becomes K due to the norm constraint, and 

e A’(i) (70) 

which hnishes the proof. 


Appendix B 

Proof of PropositionEI 

First, we show that if 

u-A/'(0,cr2), 


(71) 
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we have 


Pr (sgn(u — v) = sgn(u + v)) 


r*+oo 


1 


rl“l 


1 


-oo v2^ 
^ ^+oo 


J —oo 

o /•+00 


'U^V JQ 


(a) 2 


^+oo 


TTO-^ Jo 
2 

= — arctan — 

TT (Jy 


l-\u\ \/27rcr2 

/■|“l _ „2 
f e ^ dvdu 

e dvdu 

W ^2 

e ^ dvdw 


(72) 

^2 

e ^ dvdu (73) 


where (a) follows by letting u = w—. 
To prove Proposition [T] note that 


•\/PHxi + n = \/PhiXi + 
•\/PHx 2 + n = —y/PhiXi 


/PH'x' + n, 
■f /PH'x' + n 


where 


H' = [h2 hg • • • h^] , 

x' = [x 2 X 3 ■■■ Xk] . 
Because of the assumptions on H and Xk, we have 

'/Ph.ixi ^ CA/'(07Vr,PI at,), 

V^H'x'+ n - CA((Ow,, ((P - 1)P + CT^)lArJ. 


(74) 

(75) 

(76) 

(77) 


(78) 

(79) 

(80) 
( 81 ) 

(82) 

(83) 


Therefore, the requirement of yj^ = y 2 can be decomposed 
into 2N,- independent equations sgn(u — v) = sgn(u + v), 
which gives 

2 Nr 


Pr(yi = y 2 ) = 


arctan 


(iT- l)P + cr2 


■ (84) 
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